arXiv:1509.06787v2 [cond-mat.mtrl-sci] 25 Feb 2016 


Symmetry-selected spin-split hybrid states in Ceo/ferromagnetic interfaces 


Dongzhe Li/ Cyrille Barreteau/’^ Seiji Leo Kawahara/ Jerome Lagoute/ Cyril 
Chacon/ Yann Girard/ Sylvie Roussel/ Vincent Repain/ and Alexander Smogunov^’j^ 

^Service de Physique de VEtat Condense (CNRS UMR 3680), IRAMIS/SPEC, 

CEA Saclay, Universite Paris-Saclay, E-91191 Gif-sur-Yvette Cedex, Erance 
^DTU NANOTECH, Technical University of Denmark, 0rsteds Plads 344, DK-2800 Kgs. Lyngby, Denmark 
^Laboratoire Materiaux et Phenomenes Quantiques, Universite Paris Diderot et CNRS, 

UMR 7162, Case courrier 7021, 75205 Paris Cedex 13, Prance 
(Dated: February 26, 2016) 

The understanding of orbital hybridization and spin-polarization at the organic-ferromagnetic interface is 
essential in the search for efficient hybrid spintronic devices. Here, using first-principles calculations, we re¬ 
port a systematic study of spin-split hybrid states of Ceo deposited on various ferromagnetic surfaces: bcc- 
Cr(OOl), bcc-Fe(OOl), bcc-Co(OOl), fcc-Co(OOl) and hcp-Co(OOOl). We show that the adsorption geometry of 
the molecule with respect to the surface crystallographic orientation of the magnetic substrate as well as the 
strength of the interaction play a crucial role in the spin-polarization of the hybrid orbitals. We find that a 
large spin-polarization in vacuum above the buckyball can only be achieved if the molecule is adsorbed upon a 
bcc-(OOl) surface by its pentagonal ring. Therefore, bcc-Cr(OOl), bcc-Fe(OOl) and bcc-Co(OOl) are the optimal 
candidates. Spin-polarized scanning tunneling spectroscopy measurements on single Ceo adsorbed on Cr(OOl) 
and Co/Pt(lll) also confirm that both the symmetry of the substrate and of the molecular conformation have a 
strong influence on the induced spin polarization. Our finding may give valuable insights for further engineering 
of spin filtering devices through single molecular orbitals. 

PACS numbers: 71.15.-m, 75.70.Cn, 75.50.Pp 


I. Introduction 

The organic spintronics has become an exciting research 
field in nanoelectronics because of its flexibility, low produc¬ 
tion costs and easy functionalization. In particular, the carbon- 
based materials are promising candidates for efficient tunnel 
barrier in spintronics devices (an organic layer sandwiched be¬ 
tween two magnetic electrodes) with a large spin-dependent 
transport due to its low intrinsic spin-orbit coupling as well as 
weak hyperfine interaction. Unexpected large magnetoresis¬ 
tances and a large spin-dependent transport length have been 
reported in spin valves using Alq/^tU^ carbon nanotubeJ^, self- 
assembled molecular wireP and CeP as an organic spacer 
layer between two magnetic layers. Recently, spin coherent 
transport in spin valves based on Ceo at room temperaturP 
was also reported. 

In this context, the study of organic-ferromagnetic hybrid 
interfaces at the molecular level plays a key role to get a bet¬ 
ter understanding of the physical mechanism involved in spin 
injection and the subsequent spin-transport in organic spin 
valveP. Due to the numerous peculiar properties of Ceo (high 
symmetry, easy production, thermal and mechanical resilience 
etc), a large effort has been dedicated to the survey of Ceo- 
ferromagnetic interfaces. For instance, the strong hybridiza¬ 
tion between the d states of magnetic surfaces and the lowest 
unoccupied molecular orbitals (LUMOs) of a Coo molecule 
leads to spin-polarized molecular states close to the Fermi 
level in C6o/Cr(OOlj2!, Ceo/NiCllljSl and C6o/Fe(001 P in- 
terfaces. 

The present work is motivated by our recent study of single 
Coo molecules adsorbed on a Cr(OOl) surface at the pentagon 
site (lowest energy configuration^. There, large tunneling 
magnetoresistances (TMR) were observed by Spin-Polarized 


Scanning Tunneling Spectroscopy (SP-STS). Performing ab 
initio calculations we have attributed the TMR effect to the 
spin-splitting of one of the three LUMO orbitals close to the 
Fermi energy due to strong hybridization with the magnetic 
substrate. This orbital, labeled m = 0, is localized at the 
pentagon center and extends significantly in the vacuum. On 
another side, the two other LUMO orbitals, labeled m = ±1, 
are strongly localized on the pentagon ring and are almost not 
detectable in experiment. 

One can question, therefore, about the generality of such 
results and under which conditions (substrate nature, crystal¬ 
lographic orientation, Ceo adsorption site etc.) they can be 
reproduced or even optimized. In this paper we will thus ad¬ 
dress two main aspects: i) the strength and details of the spin- 
polarization of LUMO orbitals due to magnetic substrate and 
ii) the detectability of such spin-polarized orbitals (if any) in 
SP-STS measurements which would give rise to the TMR ef¬ 
fect. 

Our manuscript is organised as follows. In Sec. we 
present the computational and experimental methods used in 
this work. In Sec. we present a systematic ab initio study 
and show SP-STS measurements on Ceo adsorbed on two dif¬ 
ferent substrates. A particular emphasis will be made on the 
spin-resolved local density of states (LDOS) in the vaccum 
above the molecule which is the relevant quantity when com¬ 
paring to STM spectra. Finally, the conclusions will be drawn 
in Sec. 113 
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bcc-Cr(OOI) 
Pentagonal 6:6 bond 

bcc-Pe(OOI) 
Pentagon 6:6 bond 

bcc-Co(OOl) 

Pentagon 

fcc-Co(OOl) 

Pentagon 

hcp-Co(OOOl) 

Pentagon 5:6 bond Hexagon 

Eb (eV) 

3.90 

3.50 

2.32 

2.71 

2.90 

2.29 

0.79 

1.35 

1.21 

t^C-Substrate (A) 

2.05 

2.08 

2.05 

2.01 

1.98 

1.95 

2.08 

1.95 

2.08 

Charge transfer (e) 

0.50 

0.36 

0.38 

0.32 

0.30 

0.03 

0.00 

0.02 

0.08 

Ms on Ceo (/^b) 

-0.45 

-0.03 

-0.60 

-0.38 

-0.45 

-0.40 

-0.31 

-0.23 

-0.33 


TABLE I: Binding energy £^b, shortest carbon-metal interatomic distance, charge transfer to the molecule and induced spin moment Ms on 
Ceo adsorbed on cubic magnetic surfaces [bcc-Cr(001 bcc-Fe(OOl), bcc-Co(OOl) and fcc-Co(OOl)], and on a hexagonal magnetic surface 
[hcp-Co(OOOl)] with different adsorption geometries. Inequivalent adsorption geometries considered in the paper are demonstrated on upper 
panels where only the lower half of Ceo and the surface layer of the substrate are shown for simplicity. 


II. Methodology 

A. Density functional theory calculations 

Spin-polarized ab initio studies were carried out us¬ 
ing the plane wave electronic structure package Quantum 
ESPRESSCplin the framework of the density functional the¬ 
ory (DET). Generalized gradient approximation in Perdew, 
Burke and Ernzerhof parametrizatiorP^ was used for the 
exchange-correlation functional within the ultra-soft pseu¬ 
dopotential formalism. Energy cut-offs of 30 Ry and 300 
Ry were employed for the wavefunctions and the charge den¬ 
sity, respectively. The interfaces were modeled by five layer 
slabs of magnetic material with a (4 x 4) in-plane periodic¬ 
ity on which one Ceo molecule was deposited. The super¬ 
cell periodicity was increased up to (5 x 5) in the case of 
fcc-Co(OOl) and hcp-Co(OOOl) to avoid interaction between 
the buckyballs. In the ionic relaxation, the Brillouin-zone has 
been discretized by using a (4 x 4 x 1) /c-points mesh and a 
Marzari Vanderbilt cold smearing parameter of 0.01 Ry. The 
two bottom layers were fixed while other three layers of sub¬ 
strate and Ceo molecule were relaxed until the atomic forces 
are less than 0.001 eV/A. In addition, a vacuum space of about 
20 A was taken to separate two neighboring slabs in the 2 : 
direction (perpendicular to the surface) in order to avoid un¬ 
physical interactions. The electronic structure of the relaxed 
structures has been studied by using a denser (6 x 6 x 1) k- 
points mesh and two additional atomic layers were added. 


B. Experiment 

The experiments have been performed in an ultra-high vac¬ 
uum setup (P<10“^^ mbar), with a Omicron low-temperature 
STM working at 4.6 K. To obtain a Co close-packed surface, 
we have deposited by e-beam evaporation a submonolayer 
Co coverage on a Pt(III) surface at room temperature. The 
Pt(III) substrate has been cleaned by Ar ions sputtering cy¬ 
cles (I kV, P = 5.10“^ mbar) at 800 K followed by fiash an¬ 
nealing to 880 K until the Auger spectrum does not show any 
carbon contaminant. This system is known to display spin 


contrast among the different Co islands, blocked either in a 
spin up or spin down out-of-plane configuratioiP^. The spin 
polarized results have been obtained using a freshly annealed 
tungsten tip coated with around 20 nm of iron. Once the 
spin contrast has been achieved, the Ceo has been deposited 
at 4.6 K. As deposited, the adsorption configurations of the 
Ceo are numerous, even though a vertex configuration seems 
to be more frequent (cf. appendix C). Such a sample makes 
difficult an analysis of the spin contrast for different molecules 
on the same geometry. We have therefore annealed the sam¬ 
ple to room temperature, what leads to single Ceo molecules 
adsorbed on a pentagonal configuration (cf appendix C).The 
Cr(OOI) sample has been cleaned by Ar ions sputtering cycles 
(2 kV, P = 5.10“^ mbar) at 850 K. The Ceo has been deposited 
at 4.6 K onto the sample and directly measured, what allows 
to observe various adsorption geometries. It is worth noting 
that an annealing at room temperature of such a sample leads 
to a favored pentagonal configuration that has already been 
studied in detail and display a spin-split LUMO leveffl 


III. Results and discussion 


A. DFT calculations 


We have investigated various Ceo/ferromagnetic interfaces, 
i.e., different chemical elements (Cr, Ee, Co) and crystalline 
structures (cubic and hexagonal) of the substrates as well as 
different molecule adsorption sites (pentagon, hexagon, 5:6 
bond, 6:6 bond), with a particular emphasis on the hybridiza¬ 
tion between molecular levels and surface spin-polarized 
states. In Table U the main results are summarized for all stud¬ 
ied configurations. As a general result, we mention the anti¬ 
ferromagnetic (with respect to the substrate) magnetization in¬ 
duced on Ceo molecules and, on another side, the decrease of 
the spin moment for the surface atoms beneath the molecule 
due to hybridization with molecular states (see Appendix A 
for more details) - these effects have been already reported 
for Ci^land eJ^M 
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FIG. 1: (Color online) Electronic structure of a Ceo molecule ab¬ 
sorbed on bcc-Fe(OOl) by a pentagon face and in 6:6 bond geometry: 
(a) upper panels shown the total DOS of the free (but distorted) Ceo 
molecule. Lower panels show the spin-resolved density of states of 
the full Ceo/Fe system projected on the three LUMO levels of the 
isolated molecule; (b) present the spin-resolved vacuum LDOS at 5 
A above the free (distorted) molecule (upper panels) and above the 
Ceo in contact with the bcc-Fe(OOl) substrate (lower panels). Spin 
up/down curves are plotted by solid/dashed lines. 


7. Ceo on Fe bcc-(001) surface geometry 

We start our discussion with the case of a Ceo molecule ad¬ 
sorbed on a Fe-bcc(OOl) surface which we expect to behave 
similarly to the Cr-bcc(OOl) since both metals have the same 
crystalline structure and almost the same lattice parameter. 
We have considered two adsorption geometries presented in 
Fig. [2 In contrast to the case of Cr, for which the pentagonal 
configuration is the most stable, the lowest energy configura¬ 
tion for Ceo/ Fe(OOl) corresponds to the molecule bound by a 
6:6 bond to the iron surface which agrees with previous DFT 
calculations.The pentagon geometry has a higher energy of 
about 0.4 eV (see Table|I|. 

In Fig. [2 (a), for both adsorption geometries, we present in 
upper panels the density of states (DOS) for the isolated (but 
distorted, as in the contact with substrate) molecule. Lower 
panels show the DOS for the full Ceo/Fe system projected 
onto LUMO orbitals of the (distorted) Ceo molecule which 
largely dominate around the Fermi energy. Note that due to 
the structural distortion, caused by the interaction with the 
substrate, there is a breaking of the icosahedral symmmetry 
and the original three-fold degenerate Ceo LUMOs will split 
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FIG. 2: (Color online) Electronic structure of a Ceo molecule ab¬ 
sorbed on bcc-Ee(OOl) by a pentagon face and in 6:6 bond geometry: 
(a) upper panels shown the total DOS of the free (but distorted) Ceo 
molecule. Lower panels show the spin-resolved density of states of 
the full Ceo/Fe system projected on the three LUMO levels of the 
isolated molecule; (b) present the spin-resolved vacuum LDOS at 5 
A above the free (distorted) molecule (upper panels) and above the 
Ceo in contact with the bcc-Fe(OOl) substrate (lower panels). Spin 
up/down curves are plotted by solid/dashed lines. 


differently for the two adsorption geometries. 

In the ideal pentagonal geometry the LUMO states can be 
labeled by an integer m that refiects their symmetry with re¬ 
spect to the five-fold rotation axis ^sing through the centers 
of two opposite pentagonal facets!^ For deposited molecule 
this symmetry axis is lost and the only remaining symmetry 
is the refiection plane perpendicular to the surface and divid¬ 
ing the molecule into two. As a consequence, the three-fold 
degenerate LUMO levels labeled by m = 0 and m = ±1 
are split into one “even” m = 0 and two (almost) degener¬ 
ate m = ±1 (one “odd” and one “even”) levels. The infiu- 
ence of the substrate can evidently not be be summed up by 
a simple lifting of the disctete levels of the molecule: they 
will also be strongly broadened (and spin-split) due to their 
hybridization with the spin-polarized substrate states. Never¬ 
theless, it is interesting to formally divide the coupling process 
in two contributions: a first one related to the distorsion of the 
molecule and another one to the electronic hybridization. The 
splitting of the LUMO levels described above is clearly seen, 
the m = 0 being lower in energy. But the hybridization pro¬ 
duces much more drastic effect on the levels of the molecule 
which are broadened over an energy range of more than 3 eV. 
The two types of LUMOs are affected rather differently but 
still bear some resemblances: all of them are negatively polar¬ 
ized, the spin down states show a resonance below the Fermi 
level (around -0.5 eV) while spin up states have a resonance 
slightly above the Fermi level. These features seem to be more 
pronounced, however, for the m = 0 states. 

For the 6:6 bridge geometry, symmetry arguments put for¬ 
ward for the pentagonal geometry do not apply anymore since 
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no symmetry is present at all. The three LUMO states, (la¬ 
beled LUMOl, LUM02, and LUM03), are fully split as seen 
from the DOS of the isolated deformed molecule. However 
the LUMOl seems to keep a dominant m = 0 character while 
m = ±1 mix strongly. As a result the PDOS of LUMOl 
resembles somewhat the one obtained with the pentagonal ge¬ 
ometry while LUM02 and LUM03 are rather different and 
show almost no exchange splitting. 

The PDOS analysis, presented above, is very instructive 
but since one of our goals is to predict which system will 
show spin-split states measurable in a tunneling regime, typ¬ 
ically probed by a spin-polarized scanning tunneling spec¬ 
troscopy (SP-STS), the so-called vacuum local density of 
states (LDOS) is a more appropriate quantity. Indeed, in a 
simplified Tersoff-Hamann approachPI, the spin-polarized dif¬ 
ferential conductance is simply related to the spin-resolved 
LDOS of the sample at the STM tip position: 

G= ocy]n^ng(i?T,^F + eV) (1) 

cr 

where + eV) are, respectively, the spin- 

dependent tip DOS (assumed to be constant in energy) and 
vacuum LDOS of the sample (Ceo molecule deposited on sur¬ 
face) calculated at the tip position Rt above the molecule and 
at the energy corresponding to applied voltage V. 

We present in Fig. (b) the vacuum LDOS calculated 
by integrating n^{RrY,E) inside a small cubic box of size 
0.4 A at 5 A above the Ceo molecule for both spin polar¬ 
izations. For isolated molecules (upper panels), it allows to 
single out molecular orbitals with the slowest decay in vac¬ 
uum which should be accessible to the STM measurements. 
It turns out that, for both adsorption configurations, only one 
of the LUMO orbitals will mainly dominate in the tunneling 
current. In the pentagon geometry, it is m = 0 orbital which 
is essentially located in the center of the pentagon as will be 
discussed later. In the 6:6 bond configuration, it is LUM03 
orbital while the detection of the LUMOl orbital (resembling 
the m = 0) is strongly unfavored due to the rotation of the 
molecule. As a consequence, the vacuum LDOS for adsorbed 
molecules (lower panels) follows very closely the correspond¬ 
ing molecular PDOS around the Fermi level. Since, in the 
pentagon geometry, the m = 0 is clearly spin-split this results 
in large spin-polarization of vacuum LDOS leading to large 
magnetoresistance which should be visible in SP-STS exper¬ 
iment, similar to what has been measured on the Cr surfac 
In contrast, for the 6:6 bridge geometry, the LUM03 orbital 
shows almost no exchange splitting so the vacuum LDOS is 
rather spin-independent which should translate into very small 
TMR values. Note also that in this bridge configuration a spe¬ 
cific feature in LDOS is seen at -1.5 eV. It can be attributed 
to the HOMO levels which in addition show a rather modest 
spin-splitting. These states are hardly visible in the case of 
the pentagonal geometry since the electronic density of the 
HOMO states is rather low at the center of the pentagon. 

One can thus conclude that even though both adsorption 
configurations possess spin-polarized LUMO-derived states 
around the Fermi energy it is not enough to get a spin- 


polarized signal in a SP-STS measurement. It is demon¬ 
strated by the bridge geometry where the slowest decaying 
state in vacuum (and thus dominating the tunneling current) 
was found to be only weakly spin-polarized. Finally, we note 
that we have also looked at the bridge adsorption site on the 
Cr surface and have not found neither any significant spin- 
polarization of the vacuum LDOS (see Appendix B, Fig. 
(d)). 


2. Ceo on various (001) surfaces in pentagonal geometry 

Since the pentagonal adsorption site on a cubic lattice 
seems to be the best configuration to generate a strong ex¬ 
change splitting of the LUMO states that can be probed in the 
tunneling regime, it is tempting to look at other possible cu¬ 
bic substrates such as Co-bcc(OOl) and Co-fcc(OOl) and com¬ 
pare them with Fe and Cr cases discussed above. For the sake 
of completeness we present in Fig. the vacuum LDOS and 
PDOS on the m = 0 LUMO - which plays a crucial role - 
for all the four cases together. As previously mentioned, the 
Cr and Fe bear some similarities. As shown in Fig. in the 
vicinity of the Fermi level the vacuum LDOS is largely dom¬ 
inated by the LUMO states of m = 0 character which are at 
the origin of two resonances: one at —0.5 eV for the down 
spin and another one slightly above the Fermi level (0.2 eV) 
for the up spin polarization. There is also another sharp fea¬ 
ture above 1 eV but that do not show any clear spin-splitting 
and comes probably from LUMO-fI states having somewhat 
similar shape as the m = 0 LUMO orbital. The case of bcc- 
Co(OOl) is also quite similar: two spin-split LUMO m = 0 
resonances are present approximately at the same energy po¬ 
sition (the spin down feature being slightly shifted towards 
lower energy: -0.8 eV). Above the Fermi energy the broad fea¬ 
ture extending towards high energy originates from LUMO-fI 
derived states. In contrast, the fcc-Co(OOl) shows a rather dif¬ 
ferent behavior. In particular, there is a low and featureless 
electronic density below the Fermi level (down to at least -3 
eV below the Fermi level) and a modest exchange splitting 
above the Fermi level. The departure from the case of the 
bcc-(OOl) substrates can probably be attributed to the dens¬ 
est atomic packing of the fcc(OOl) surface: the interatomic 
distance of the 2D square lattice of fcc-Co(OOl) (a = 2.49 
A) is much smaller than the one of bcc-Co(OOl) (a = 2.84 
A). Therefore, more carbon atoms of the Ceo molecule are in¬ 
volved in a bonding with the substrate but the strength of the 
hybridization due to a given metal-carbon bond is weaker. It 
results in an overall binding energy similar to the case of the 
bcc-Co(OOl) surface (see Table I) but since the vacuum LDOS 
is still largely dominated by the contribution of the m = 0 
molecular orbital (which is less hybridized) a smaller spin¬ 
splitting is found. We can summarize that the spin-splitting of 
the m = 0 LUMO becomes less and less pronounced when 
passing from the bcc-Cr to fcc-Co accompanied by a contin¬ 
uous disappearance of related spin-split peaks in the vacuum 
LDOS. 

To further rationalize this behavior we have tried to quan¬ 
tify the molecule distortion and relate it to the LUMO splitting 
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FIG. 3: Electronic structure of a Ceo molecule adsorbed on top of cubic magnetic surfaces in pentagonal ring geometry: (a) bcc-Cr(001 (b) 
bcc-Fe(OOl), (c) bcc-Co(OOl) and (d) fcc-Co(OOl). Upper panels: Spin-resolved vacuum local density of states (LDOS) at 5 A above the Ceo 
molecule. Lower panels: Density of states of the Ceo/substrate projected onto the molecular orbital LUMO m = 0 of the deformed molecule. 
Spin up/down curves are plotted by solid/dashed lines. 


in the four cases discussed above (Fig. |^. The deformation of 
the molecule in contact with the substrate is evaluated by com¬ 
parison with the free molecule in a direction perpendicular and 
parallel to the surface, namely, ALz and ALx, respectively. 
As shown in Fig. |^(a), the flattening of the molecule (ALz) is 
much more important than its lateral expansion (ALx) which 
can be neglected. Interestingly, the deformation is basically 
proportional to the binding energy (see Table for the bcc 
surfaces while almost no distortion is found for the fcc-Co sur¬ 
face although the interaction is substantial. Once again this is 
probably due to the more dense packing of the surface atomic 
layer of the fee structure. 

The level splitting of three-fold degenerate LUMO levels 
goes in line with the amplitude of the deformation as demon¬ 
strated by Fig. 1^ (b). The main out-of-plane deformation 
ALz, preserving the flve-fold rotational axis, will lift the de¬ 
generacy between m = 0 and m = ±1 orbitals, with the gap 
AE roughly proportional to ALz • The m = 0 has a circular 
shape extending out of the molecule at the center of the pen¬ 
tagon while m = ±1 states have a strong localization on the 
pentagon ring. The additional small (and slightly anisotropic, 
ALx ~ ALy) in-plane distortion, and the electronic coupling 
to the cubic substrate will break the flve-fold rotational sym¬ 
metry leaving only mirror symmetry with respect to the plane 
yz (see Fig. |^(b)) passing through the molecule and perpen¬ 
dicular to the surface. The m = 0 is of “even” symmetry 
while m = ±1 states will be slightly split into one state of 
“odd” and one of “even” symmetry. The important point is 
that this latter “even” state can admix to the m = 0 since they 
have the same symmetry, destroying its circular shape, and the 
amount of this admixture depends on the gap AE - the bigger 
AE the smaller the admixture. 

Since the m = 0 state is crucial in producing the spin-split 
features in vacuum LDOS, as has been shown above, it is cru¬ 
cial to keep it as pure as possible. Therefore, the larger ALz 


the more pronounced are the two spin-split peaks in the vac¬ 
uum LDOS which explains the general tendency observed in 
Fig-lH 


3. Ceo on hexagonal magnetic surface 

Finally let us consider another important substrate with a 
different symmetry: hcp-Co(OOOl). The most stable relaxed 
atomic conflguration of the Ceo was found to be the 5:6 bond 
(see Table I) where the molecule is bound to a surface Co 
atom by a pentagon-hexagon edge. Nevertheless since it is 
rather easy for this system to be trapped in a local metastable 
minimum we have considered the two other conflgurations, 
namely pentagonal and hexagonal. In Fig. we plot the 
vacuum LDOS of Ceo absorbed on hcp-Co(OOOl) for these 
three different adsorption sites. As a general result, in all 
three cases the vacuum LDOS is very low around the Fermi 
level and show almost no spin-polarization. Interestingly, in 
the pentagonal geometry the PDOS analysis (see Supplemen¬ 
tary materials) reveals again slight spin-polarization of the 
LUMOl orbital around —0.5 eV. However, due to symme¬ 
try mismatch with the substrate, ah the LUMOs are mixed 
and decay almost equally to the vacuum which results in very 
weak spin-polarization of LDOS above the Ceo molecule. In 
the hexagonal geometry, on the contrary, only LUMOS or¬ 
bital has a large extension in the vacuum, but it is not spin- 
polarized (see Supplementary materials) which results again 
in weak spin-polarization of the vacuum LDOS. We can con¬ 
clude, therefore, that for observing a low-bias TMR signal in 
Ceo/ferromagnet systems it is crucial to have only one low- 
decaying LUMO state which, in addition, should be spin- 
polarized due to hybridization with ferromagnetic substrate. 
This was found to happen only for m = 0 LUMO orbital 
in the pentagon adsorption geometry on cubic magnetic sub- 
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FIG. 4: (Color online) (a) Ceo deformation due to the interaction 
with a (001) surface of cubic materials. AL^ (AL^) quantifies 
the deformation of the adsorbed molecule (with respect to the free 
molecule) in the direction z perpendicular (x parallel) to the surface. 
A schematic illustration of the deformation is also presented on the 
top panel, (b) Energy level diagram of LUMO levels for the free but 
distorted Ceo molecule (after ionic relaxation of the full system). The 
3 -fold degeneracy of the Ceo LUMOs is lifted by the molecule defor¬ 
mation into a single m = 0 state and two almost degenerate m = ±1 
states which can be classified as ’’even” or ’’odd” with respect to the 
mirror plane yz. The spatial representation of the three LUMO levels 
is also presented, the isosurfaces of positive and negative values are 
shown in red and blue. 


strates. 

It should be finally emphasized that if the contact regime 
is concerned (rather than the tunneling regime), which could 
be realized in spin valve devices using a “sandwich” trilayer 
geometry, then the vacuum decay rate of different LUMO 
orbitals is no longer of importance and all the LUMOs are 
expected to contribute to the current. In this case the stan¬ 
dard spin-polarized PDOS analysis presented in Fig.[2(a) for 
Ceo/Fe system is already sufficient. In particular, it suggests 
that in both Ceo/Fe adsorption geometries a spinpolarized 
LUMO orbital shows up, so that the spin-polarized current 
should be observed irrespective of the details in Ceo/Fe inter- 
facethe result recently reported for the Fe-bcc(001)/C60/Fe- 
bcc(OOl) junction^. 


hcp-Co(OOOl) 
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FIG. 5: Calculated spin-resolved vacuum LDOS at 5 A above the 
Ceo molecule deposited on hcp-Co(OOl) surface with various adsorp¬ 
tion sites: (a) a pentagon ring , (b) a 5:6 bond, (c) a hexagon ring. 
Spin up/down curves are shown by solid/dashed lines. 


B. SP-STS measurements 


In order to confirm that the spin polarization probed by 
SP-STM above a Ceo molecule is very sensitive to the sub¬ 
strate crystal structure and to the adsorption geometry of 
the molecule, we have measured single Ceo adsorbed on 
Co/Pt(l 11) and on Cr(OOl). It is worth noting that it is gener¬ 
ally difficult to achieve a quantitative comparison between ex¬ 
perimental and theoretical results for several reasons. Firstly, 
the experimental systems can display complex structures (dis- 
commensuration lines on Co/Pt(lll), cf. appendix C) and lo¬ 
cal defects that cannot be taken into account in calculations as 
they would require very large unit cells and a too long compu¬ 
tational time. Secondly, the LDOS calculation in vacuum on 
a single point (Tersoff-Hamann model) does not give always 
an accurate description of STM results. Molecular orbitals are 
generally broader in experiments than calculated and at differ¬ 
ent energieJJ^. It should be even worth for SP-STM where the 
tip polarization is likely to play an important role in the quan¬ 
titative analysis of spin polarized spectra. Fig. (a) shows 
the results for Ceo on Co/Pt(lll). The image in inset shows 
several one monolayer high Co islands with a color code dis¬ 
playing their differential conductance at -1 V which is known 
to display a spin contrasP^. The spectra averaged over two 
opposite spin direction islands (yellow and green) show no 
feature above Ep and broad peaks in the negative range, from 
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0 to -1.3 V, that have been ascribed to electronic surface res¬ 
onances of CcP. The blue and red spectra have been taken 
above six single Ceo molecules, all adsorbed in a pentagonal 
configuration, as checked by high resolution images, but in 
contact either with a spin down Co island (blue spectra) or 
with a spin up Co island (red spectra). All these six spec¬ 
tra are very similar and typical of Ceo adsorbed on a close- 
packed surfacJi^ with a peak associated to the HOMO level 
at -2 V, to the LUMO at 1.1 V and to the LUMO-f 1 at 2 V. 
The important result is that we cannot observe any significant 
difference between the molecules adsorbed on spin down and 
spin up islands. The comparison with theoretical calculations 
can only be qualitative as the real structure of the Co islands 
is very complex with different stacking area, discommensu- 
ration lines and surface dislocation^^. However, the absence 
of measurable spin polarization is in good agreement with the 
conclusion of Sec. |III A^ where all the adsorption geometries 
on Co(OOOl) show a very weak spin polarization. 

Fig.|^(b) focus on the case of Ceo molecules adsorbed on 
a roughly hexagonal configuration on a Cr(OOl) surface. A 
closer look at high resolution images, such as images in inset, 
shows that the studied molecules are in between an hexag¬ 
onal and a 6:6 adsorption geometry. The typical spectra of 
the raw Cr surface for spin down (yellow) and spin up ter¬ 
races (green) show a sharp features close to the Fermi level 
that has been associated to ^ s — pz — dz^ surface state and 
an electronic state around 1 V, both displaying spin contrast 
on the conductance image^SEHl^ Xwo spectra associated with 
a molecule lying on a spin down terrace (up images in inset) 
and a molecule on a spin up terrace (down images in inset) 
are shown respectively in blue and red. These spectra are very 
different from the ones measured on molecules adsorbed in 
a pentagonal configuratioJ^. More specifically, they do not 
show a spin-split LUMO orbital, in good agreement with the 
theoretical findings. Mainly two LUMO states at 0.7 V and 
1.1 V are observed, with only a slight change in intensity for 
spin up and spin down configurations for the LUMO-i-1. This 
can be compared with the results of Fig.[^(d) (Appendix B), 
having in mind that the experimental adsorption geometry is 
not exactly the same as the calculated one. The quantitative 
agreement is pretty good with a LUMO level insensitive to the 
spin polarization of the substrate and a LUMO-fI level that is 
slightly different for the spin up and spin down configurations. 

Finally, those experimental measurements on Co(OOOl) 
and Cr(OOl) confirm that the vacuum LDOS above the Ceo 
molecules only shows a clear spin polarization and spin-split 
orbital in a pentagonal adsorption geometry on the (001) sur¬ 
face of a bcc. Other adsorption geometries or different sub¬ 
strates like Co(OOOl) are less favorable to observe such spin- 
split molecular orbitals. 


IV. Conclusions 

To conclude, we have investigated systematically from first- 
principles, the spin-polarized hybrid states of Ceo deposited 
on ferromagnetic surfaces such as bcc-Cr(OOl), bcc-Fe(OOl), 
bcc-Co(OOl), fcc-Co(OOl) and hcp-Co(OOOl). As a general 



—I-■-1-■-1-■-1-■-r 

- 2-1012 

bias voltage (V) 



FIG. 6: (Color online) (a) Differential conductance as a function of 
the tunneling voltage recorded on Co/Pt(lll) and on Ceo adsorbed 
on a pentagonal configuration. Green and orange spectra are the ref¬ 
erence taken on up spin and down spin Co islands. Blue and red 
spectra are taken on top of six different Ceo molecules adsorbed re¬ 
spectively on up spin and down spin Co islands. In inset, a 3D to¬ 
pographic STM image (50 x 50 nm^) of Ceo molecules adsorbed on 
Co islands on Pt(lll). The color code is the conductance at -1 V, 
showing the spin contrast in the Co islands, (b) Differential conduc¬ 
tance as function of the tunneling voltage recorded on Cr(OOl) and 
on Ceo adsorbed on an hexagonal configuration. Orange and green 
spectra display the signal of two adjacent Cr terraces, i.e. down spin 
and up spin. The blue and red curves show spectra averaged over 
a whole Ceo molecule adsorbed on down spin and up spin terrace 
respectively. For clarity, those latter curves have been magnified by 
a factor three. In inset, conductance images of the corresponding 
Ceo molecules, taken at the energies of the molecular levels (0.7 and 
1.1 V). 


feature, a strong chemisorption of the buckyball on ferromag¬ 
netic surfaces leads to a remarkable drop of the spin moment 
for magnetic surface atoms and an induced negative spin mo¬ 
ment for the Ceo molecule. It was found that the degree of 
spin-polarization of the Ceo LUMO depends strongly on both 
the adsorption geometry and the symmetry of the surface. Due 
to symmetry matching between the Ceo LUMO states with 
the underlying surface and a non-negligible distortion in per¬ 
pendicular direction to the surface, the 3-fold degeneracy of 
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LUMOs of free Ceo could be lifted energetically and spa¬ 
tially. As a result, a large spin-polarization of a single hy¬ 
brid orbital is only achieved if a Ceo molecule is adsorbed 
by a pentagon face on cubic surfaces, such as bcc-Cr(OOl), 
bcc-Fe(OOl) and bcc-Co(OOl), and is related to the spin-split 
m = 0 LUMO which is strongly localized at the centre of the 
pentagon ring. In contrast, the adsorption on the hexagonal 
hcp-Co(OOOl) surface leads to very small vacuum molecular 
induced spin polarization. Our theoretical results are qualita¬ 
tively confirmed by SP-STS measurements that show no mea¬ 
surable spin polarization for Ceo on Co/Pt(lll) and no spin- 
split LUMO orbital for a non pentagonal adsorption geometry 
on Cr(OOl). Understanding the mechanism of spin-polarized 
hybrid states at the interface is an essential ingredient for the 
engineering of spin filtering in carbon-based spintronics de¬ 
vices and we expect that this work will help for their future 
rationalized design. 
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Appendix A: Magnetic spin moment 

The calculated magnetic spin moment of C6o/bcc-Fe(001) 
and C6o/hcp-Co(0001) is shown in Fig. |^for the lowest en¬ 
ergy configurations, a strong modification for both surface and 
adsorbate have been found with the concomitant interaction 
bewteen molecule and substrate. 

We present the spin moment of the first layer of surface 
(see the left panel of Fig. [^, in both Fe and Co surfaces we 
notice the signaficant reducement of the spin moment in the 
vicinity of molecule. This drop of spin moment in surfaces 
originates from the hybridization between d states (particu- 
lally pronounced for the out-of-plane extended d orbitals) of 
surfaces and 7r-molecular orbitals. In particular, the largest 
decrease of spin moment occurs for surface atom just be¬ 
low molecule (marked in red) of Mg = 1.83 and 1.26 /iB for 
C6o/Fe(001) and C6o/Co(0001), respectively, which means a 
drop of about 40 % and 30 % compared to the clean Fe(OOl) 
and Co(OOOl) surfaces. Additionally, the spin moment of this 
particular atom becomes even smaller than the bulk value of 
2.30 and 1.64 /iB for bcc-Fe and hcp-Co, respectively. 

In addition, the induced magnetic moment of Ceo is polar¬ 
ized negatively and mainly localized around to the interface 
(see the right panel of Fig. [7]). The spin moment appears 
to oscillate when moving away from the interface and con¬ 
verge to the expected value zero finally. We estimated the total 
magnetic moment of Ceo by summing over the local spin mo¬ 
ments calculated by projecting the spin-resolved Kohn-Sham 
states onto the atomic wavefunctions, it was found to be of 
about -0.38 and -0.23 /iB for C6o/Fe(001) and C6o/Co(0001), 
respectively. This transferred spin moment in the Ceo origi¬ 


nates from the three electronic levels of LUMO close to the 
Fermi level. 

Appendix B: Spin-polarized vacuum LDOS 

We present in Fig. the spin-polarized DOS of a Ceo 
molecule absorbed on hcp-Co(OOOl) with three different ab¬ 
sorption sites and on bcc-Cr(OOl) with 6:6 bond. Although the 
three LUMO orbitals are generally polarized we do not find 
any significant spin-polarization in the vacuum LDOS around 
the Fermi energy in all the cases. 

Appendix C: Adsorption geometry of Ceo on Co/Pt(lll) 

We present experimental details on the geometry of adsorp¬ 
tion of single Ceo molecules on Co/Pt(lll), as deposited on 
the sample at a temperature of 4.6 K [Fig.[^(a) and (b)] and 
after an annealing of the sample at room temperature during 
few minutes [Fig.[^(c) and (d)]. As can be seen in Fig.[^(a) 
and (c), Co islands are quasi triangular and single layer high. 
One can observe brighter lines in between the Ceo molecules 
that have been ascribed to discommensuration lineJ^^, sepa¬ 
rating different stacking area (Co in hep or fee stacking sites 
with respect to the Pt(lll) surface layer). This makes the ex¬ 
perimental system rather complex and prevents to ascribe a 
precise adsorption site for the Ceo molecules. Nevertheless, 
conductance images recorded at the energies of the different 
molecular states reveal unambiguously their adsorption ge¬ 
ometry. Fig. (b) is a conductance image recorded at 1 V 
that shows the local symetry of the different molecules. For 




FIG. 7: (Color online) Magnetic spin moment of (a) C6o/bcc-Fe(001) 
6:6 bond, (b) C6o/hcp-Co(0001) 5:6 bond. Left panels: spin moment 
in /iB of the surface atoms, the surface atom just below the bond 
is marked in red. The magnetic moment of bulk and clean surface 
are also indicated. Right panels: The distribution of induced spin 
moment in the Ceo with respect to the number of atoms, note that 
trajectory for numbering the atoms beginning from the carbon atoms 
at the interface and moving away from the interface. 
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(a) hcp-Co(OOOl) pentagon 


(b) hcp-Co(OOOl) 5:6 bond (c) hcp-Co(OOOl) hexagon 


(d) bcc-Cr(OOl) 6:6 bond 



FIG. 8: (Color online) Electronic structure of a Ceo molecule absorbed on top of hcp-Co(OOOl) with (a) a pentagonal ring, (b) 5:6 bond, (c) a 
hexagon ring geometry and (d) on top of bcc-Cr(OOl) with 6:6 bond. On upper panels the total DOS and the spin-resolved vaccum LDOS of 
the free (but distorted) Ceo molecules are shown. On lower panels the spin-resolved vaccum LDOS and PDOS on three LUMOs of the full 
Ceo/ferromagnet system are presented. Spin up (down) data are plotted by solid (dashed) lines. Spin polarization on each LUMO level is also 
indicated. 


example, three of them are nearly identical in a vertex con¬ 
figuration. After an annealing at room temperature, Fig. 
(c) shows that the molecules are surprisingly still isolated. A 
close look at the associated conductance image showing the 
LUMO state [Fig. Fig.|^(d)] shows that the adsorption geom¬ 
etry has changed to pentagonal. Indeed, the dot and ring like 


structure observed in Fig.[^(d) is very close the LUMO level 
observed and calculated for a pentagonal adsorption geometry 
on Cr(001 j3. It is worth noting that the favored pentagonal ad¬ 
sorption geometry is very unusual on a close-packed surface 
and is probably due to the complex Co structure on Pt(l 11). 
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FIG. 9: (Color online) (18 x 18) nm^ STM images of single CeO 
molecules (white round shapes) on Co/Pt(l 11) (a) topography image, 
as deposited at 4.6 K. (b) associated conductance image at the tunnel¬ 
ing voltage 1 V corresponding to the experimental LUMO level, (c) 
topography image recorded at 4.6 K after an annealing of the sample 
at room temperature, (d) associated conductance image at the tun¬ 
neling voltage 1 V, close to the LUMO level, (c) and (d) have been 
recorded with a FeAV spin polarized tip, on the same molecules that 
are measured in Fig.|^(a). 





